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^'. Abstract: 

Oh! 

(-H I In an effective Lagrangian approach to QCD we nonperturbatively calculate an analytic 

approximation to the decay rate of a false vacuum per unit volume, T/V. We do so for 



both zero and high temperature theories. This result is important for the study of the 
early universe at around the time of the QCD phase transition. It is also important in 
_c^ I order to determine the possibility of observing this false vacuum decay at the Relativistic 

Heavy Ion Collider (RHIC). Previously described dramatic signatures of the decay of false 
vacuum bubbles would occur in our case as well. 



1 Introduction 

Effective Lagrangian techniques have proven to be very valuable in Quantum Field Theory. 
There are two main types of effective Lagrangian formulations currently in use. The first 
type is the Wilsonian effective action which describes the low energy dynamics of the 
lightest particles in the theory. The second type of effective Lagrangian is defined as 
the Legendre transform of the generating function of connected Green functions. This 
formulation implements at the Lagrangian level certain anomalous Ward identities relating 
vacuum condensates of the fields and has been referred to as the anomalous effective 
Lagrangian [^ . This type of approach is very useful for studying vacuum properties and 
was first applied to supersymmetric (SUSY) gauge theories for both Yang-Mills theory 
0] and full super-QCDiP and more recently to their non-supersymmetric versions 0] |p. 
The second type of effective Lagrangian approach is the one used in this paper. 

An effective Lagrangian, or more precisely the potential part, for non-supersymmetric 
Yang- Mills theory was described in 0] and for non-supersymmetric QCD in P). This 
effective potential for QCD for Nf light flavors and Nc colors was described in [Q and 
analyzed in detail in |jl|]. This study is interesting for several reasons. First, it provides a 
generalization of the large Nc Di Vecchia-Veneziano-Witten Effective Chiral Lagrangian 
(ECL)Q Q for arbitrary Nc after integrating out the heavy "glueball" fields. Further- 
more this approach to the derivation of the WW ECL fixes all dimensional parameters 
in terms of the experimentally measurable quark and gluon condensates. Second, this 
effective Lagrangian approach makes it possible to address the problem of 6'-dependence 
in QCD. The problem of 6'-dependence is also directly related to the problem of a realistic 
axion potential since the axion comes from giving the 6'-parameter the status of a dy- 
namical field and the axion potential comes from the ^-dependence of the vacuum energy, 
Evac{(^)- Third, metastable vacua in the effective potential may provide a mechanism for 
baryogenesis at the QCD scale |]^. Finally, the previously mentioned metastable vacua 
have dramatic signatures and would potentially be observable at RHIC. The final two 
reasons provided the original motivation for the present work. 

It is important to stress that the Anomalous Effective Potential is simply a candi- 
date for an effective potential for QCD. We do not wish to imply that effective chiral 
Lagrangians are in any way wrong. In particular, we are definitely not saying that the 
Di Vecchia-Veneziano-Witten Effective Chiral Lagrangian is wrong. In fact one of the 
reasons for studying this theory is that it agrees with the WW ECL in the large Nc 
limit. This effective Lagrangian is a generalization of the WW to the case of finite Nc. 

In the detailed study of the effective potential for QCD [|T| it was found that the ef- 
fective potential for the phases of the chiral condensate exhibits cusp singularities as a 
result of topological charge quantization^. These cusp singularities act as potential barri- 
ers separating metastable vacuum states from the true physical vacuum^. The existence 
of these metastable vacua leads to the well known phenomenon of false vacuum decayf^ 
which may have consequences in axion physics [^ |]ll, baryogenesis |jl|] and/or many other 

^This effective Lagrangian is of the second type described above. 

^For specific values of parameters of the theory which cannot be definitely determined from our present 
knowledge 

■^For certain parameter values these metastable vacua occur at every value of 6. 



unexplained phenomenon which may occur during/after the QCD phase transition. Do- 
main wall solutions interpolating between a metastable vacuum and the true vacuum were 
presented in |]T]] which will be used in this paper to calculate the decay rate of the false 
vacuum. 

These metastable vacua are somewhat controversial so some comment is required. 
First we would like to point out that nontrivial vacua have been shown to exist in Yang 



Mills theories in the large Nc limit using the AdS/CFT correspondence [T^] and the same 



phenomenon was observed[ll| in the analysis of soft breaking of supersymmetric models. 



Based on the original analysis of the WW ECL[y] and recent work||12|| it is tempting to 
conclude that metastable vacua only exist in QCD for 6 ^ it. This would contradict our 
claim that there are a series of metastable vacuaQ for ^ = and certain values of some 
integer parameters of the theory. However, as was stated above, the effective Lagrangian 
we use is a generalization of the WW ECL to finite Nc and it is not unreasonable to 
believe that this generalization can lead to new features. The effective potential we use 



also represents a generalization of the effective potential of [T^ because of the inclusion 
of the f]' and the inclusion of two integer parameters, p and q, about which there is some 
controversy. The first comment we would like to make is that if we integrate out the 
singlet rj' field from our form of the effective potential we obtain the effective Lagrangian 
of [|T2|. The second comment that we would like to make is that if we choose q = 1 and 
p = Nc, as some approaches to determining their values would suggest, then there are no 
metastable vacua at 6* = and all our results would agree with [|r^. However, the values 
q = 8 and p = llNc — 2Nf arise in a number of different approaches [|l|. In this case 
we observe a series of metastable vacua and cusps in the potentials for the chiral fields 
even at ^ = 0. In any case it appears that the existence of nontrivial vacua is a general 
phenomenon for gauge theories in the strong coupling regime. The main goal of this paper 
is to develop techniques which could be useful in the study of this type of phenomenon. 
The decay rate of false vacua in large Nc Yang Mills theory was estimated in |l^ . 



Our calculation of the decay rate of a false vacuum differs from |T^ because it is valid for 
finite Nc andecause the heavy glueball degrees of freedom have been integrated out in our 
approach, while their calculation derives entirely from gluodynamics. For more comments 
on the elimination of heavy degrees of freedom see below. The decay of a false vacuum was 
estimated in |T^ in the case that 9 ^ n and Nc is finite. Our calculation differs from []T^ in 



the inclusion of the singlet rj' field and choice of parameters g = 8 and p = llNc~2Nf. As 
well, both of these estimates only use the semiclassical approximation, while we determine 
the effect of the first quantum corrections. The form of the semiclassical decay rate in 



our case is identical to that of |T3|, [T^. We make no numerical comparisons as the false 
vacua involved are all different. 

These false vacua and domain walls could lead to many interesting consequences in 
the evolution of the early universe at around the time of the QCD phase transition. One 
example of this is related to baryogenesis and dark matter, and is described in [0. The 
zero temperature decay rate of the false vacua calculated in this paper is relevant to this 
particular application. As well, these metastable states can hopefully be experimentally 
studied at RHIC and the high temperature decay rate calculated in this paper would be 



*The number of these metastable vacua goes to infinity as Nc goes to infinity. 



relevant to this research. Bubbles of this false vacuum would display CP odd signatures 



such as those described in |]IJ] where the large N^ limit was assumed. 

The purpose of the present paper is to determine how the decay rate per unit volume, 
T/V, of the false vacuum depends on physical parameters. We will use several approx- 
imations along the way to determine what we believe to be the dominant contribution 
within a factor of about the order of unity. It should be noted that while this is an ap- 
proximation it is nonperturbative in the sense that it should contain contributions from 
all perturbative diagrams. 

Fermions (nucleons) could drastically change the results but since they also drastically 
increase the difficulty of the calculation we will leave them out in our first approximation. 
The decay rate of the metastable state could only be decreased by the inclusion of fermions 
and thus our calculation is an estimate of the upper bound on the decay rate. 

As well there is the consideration that intrinsic heavy degrees of freedom might play a 



role. It has been suggested |T^ that the cusps in the potential are due to the integration 
out of heavy degrees of freedom and that the presence of cusps invalidates the construction 
of the domain wall from the effective Lagrangian. This would mean that the heavy degrees 
of freedom must be included for the complete calculation. However, as above, the decay 
rate can only be decreased by the inclusion of these heavy degrees of freedom and we 
ignore their effects in our estimate of the upper bound on the decay rate. 

Even without consideration of these interesting applications, our method of calculation 
of the determinantal prefactor is useful as an alternative to previous methods. Previous 



calculations [|T6[ ^T^ [|T8| [|T9| of bubble nucleation rates use a particular method for 



calculating the determinant ratio of operators of the form: 

M = -V^ + uj^ + aV{r), (1) 



involving a theorem from |jT6[. We prefer to use a more direct approach. Our procedure 
provides a method for obtaining both analytical approximations and exact numerical 
calculations for this determinantal prefactor. Our numerical approximation involves only 
numerical integration unlike the method of [^ which also involves numerical solution 
of differential equations and might not be reliable in some cases involving non-smooth 
perturbation potentials. Our methods are more along the lines of [^ but are sufficiently 
different to constitute independent results. As well we have used the same methods to 
calculate the decay rate in the zero temperature theory which has not been done before. 
The paper is organized as follows. In Sect. 2 we review the structure of the effective 
potential and the domain wall solution|jl|. In Sect. 3 we determine the semiclassical ap- 
proximation to the decay rate for both zero and high temperature. We evaluate the first 
quantum corrections at zero temperature in Sect. 4 and at high temperature in Sect. 5. Fi- 
nally in Sect. 6 we discuss the implications of our results for baryogenesis and dark matter 
and for observations of parity odd bubbles at RHIC 



2 The Effective Potential 

The effective potential |E[|l| for QCD is: 
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where the lightQ matter fields are described by the unitary matrix, Uij, corresponding to 
the phases of the chiral condensate: 
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where A" are the Gell-Mann matrices of SU{Nf), vr" is the octet of pseudoscalar fields 
(pions, kaons and the eta meson) and r]' is the SU{Nf) singlet pseudoscalar field. M = 
diag{mi\{^^'$^)\) , V is the 4-volume and the integers p and q are relatively prime param- 
eters. The values of the parameters p and q are not known as different proposals for their 
determination lead to different values |Q. We will not use specific values of p and q aside 
from the restriction that q y^ 1. We would like to mention that the values p = llNc — 2Nf 
and q = 8 arise in a number of different approaches and that q/p ~ 1/Nc for the U(l) 
problem to be solved. E = {bas/ {3271)0"^) where b = ^^c - |A^/ is the Gell-Mann - 
Low /5-function of QCD. The physical input to this equation are the values of the vacuum 
condensates, (a^/vr G^) = 0.012 GeV^ and (^^) = (240 MeV)^ and the quark masses. 
We use the values /^ = 132 MeV and /^/ = 86 MeV. 

We will take this potential as our starting point motivated by four of its most important 
properties (for details see [|I|): 

i) it correctly reproduces the WW effective Chiral Lagrangian in the large Nc 
limit. 

ii)it reproduces the anomalous conformal and chiral Ward identities of QCD. 

iii)it reproduces the known ^-dependence for small 6'-angles [^] but leads to 27r peri- 
odicity in 6 of physical observables. 

iv) the related effective Lagrangian for pure gluodynamicsQ] has the nice property: 
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(5) 



which was advocated earlier by Veneziano for the f/(l) problem to be resolved p2 



This effective potential is not representable by a single analytic function in the \/^ -^ oo 
limit. The thermodynamic limit selects, at each value of 6' — ilogdetf/, one particular 



^Note that the -q' is not really very light, but it enters the theory in this way. 
^Note that mixing of the flavor eigenstates is ignored at this level. 



branch {ie. a particular value of the integer /) and cusp singularities occur where the 
branches coincide. 

The effective potential for the chiral phases of the matter fields, becomes: 
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(7) 



when we take, U = diag{expi(f)g) []. This is now a piecewise smooth potential for the 
phases of the chiral condensate with cusp singularities. 

At this point some comment is required about the removal of the singlet rj' field from 
the effective potential. If we minimize the potential for / = with respect to the singlet 
field, J24>i^ we clearly obtain the solution: 



E' 



9. 



Integrating out the singlet field amounts to substituting this solution in the effective 
potential leading to: 
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(9) 



which upon the substitutions, 0i -^ 6/3 — a, 02 -^ ^/3 — f3 and Mi -^ m becomes exactly 
the potential given in Equation (2) of [|12|. Therefore, as we asserted in the Introduction, 
integrating the singlet field out of our effective potential we obtain the same effective 



potential as in [^ . 

For g 7^ 1 there exist metastable vacuum states in addition to the lowest energy 
physical vacuum which leads to false vacuum decay. For our purposes we will use the 
simplified setting where ^ = and Nf = 3, equal quark masses rrii = A MeV 0, and equal 
chiral phases (pi = cj). This amounts to studying only radial motion in the 0-space. We 
analyze the problem in the spirit of Ref.p] and only consider transitions between the 
lowest energy metastable state and the physical vacuum. The results should be easily 
generalizable to other transitions. 

For ease of calculation we rescale and shift the chiral field -^ {2 / f^, J N f)(l) — n / {qN f) 
in order to have the standard normalization of the kinetic term and a symmetrized form 
of the potential. The effective potential for ^ = becomes 
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^This is not a restriction since the quark mass matrix can always be diagonalized. 
®Note that for q ^ I there is no phenomenological sensitivity to the values of the light quark masses 
and we are allowed to choose the quark masses to be equal. For details see m. 
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Figure 1: Effective potential at ^ = for equal chiral phases 



The effective potential (|TOD has a global minimum at 0+ = {TxfT^)/{2qJNf) and a local 
minimum at 0_ = —{'n-fT,)/{2qJNf), with a cusp singularity between them, (see Fig. |l]). 

The minima are interpreted as two vacua separated by a high potential barrier (~ G^) 
which is fairly wide, while the energy splitting, AE between the states is fairly small in 
comparison. 
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Therefore we can use the thin wall approximation [Q in our calculations. The domain 
wall solution in this approximation corresponding to the effective potential (p!0[) is: 
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where Xq is the position of the center of the domain wall and 
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is the inverse width of the wall. The solution (0) is shown as a function oi x — Xq in Fig. 
(0). Its first derivative is continuous at x = Xq, but the second derivative exhibits a jump. 
The wall surface tension is given by: 
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Figure 2: Domain wall profile. 



In what follows we use this domain wall solution in our calculation of the decay rate 
per unit volume of the false vacuum for both zero and high temperature field theory. The 
next section describes the semiclassical calculation. 



3 Semiclassical Theory 

The fate of the false vacuum was discussed in [§ and the expression of T/V has the form: 

r/r = Ae-^/^[i + o(;i)]. (15) 

The semiclassical approximation at zero temperature tells us that B is given by the 
Euclidean action of 6: 



B = Sa 
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i(|y^4w)^ + w) 



(16) 



where U{(j)) is the potential, W^jf, for the chiral phases described in the previous section 
neglecting the energy difference, Ai5, between the two vacua: 
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In order for this to be finite we must have limj.^oo0(^) = 0- = T^f-nl'^'iJNf where 



r = \Jx'^ + y^ + z^ + t'^. The solution of this problem is the four dimensional equivalent 
to the solution which Coleman calls "the bounce" [|[. It describes a bubble of the true 



vacuum in the false vacuum at the origin which forms, grows to a maximum size and then 
shrinks to nothing again (^e. an 0(4) invariant bubble). In the thin wall approximation 
the bounce solution is: 
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and the value of B is calculated to obtain[Q: 
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where a is the domain wall energy and: 
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is the Euclidean action of the 4D bounce solution. 

The thin wall approximation is valid because the radius of the 4D bubble, which is 
found to be i? = 3a/ AE by minimizing the value 5*4 [0;,] , is much larger than the width 
of the domain wall, l//i. 

For finite temperature QCD the semiclassical approximation is slightly different. At 
sufficiently high temperature the bubble solution becomes a stable 0(3) invariant bubble 
with radius p3[|: 

RiT) = ^. (21) 

where the temperature dependence of the domain wall energy, Si{T), is not known. In 
this case the calculation of B gives EH]: 
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It should be noted that this decay rate is for the ground state of the metastable well. The 
decay rate for excited energy states above the metastable vacuum via thermally activated 
transitions while similar in form is not the same[^. 

This completes the semiclassical analysis of the decay rate. In the next section we 
calculate the quantum corrections to the decay rate in the zero temperature theory. 



4 Quantum Corrections at Zero Temperature 

The quantum corrections at zero temperature correspond to the coefficient A in Eq. (p!5|) p5[] : 
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The spectrum of the operator in the numerator consists of a discrete spectrum with 
zero and negative eigenvalues and a continuous positive eigenvalue spectrum starting at 
^2 ^ [/"{(p-). These two parts of the spectrum must be analyzed separately and it can 
be shown that this factor separates into three parts: 
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The first term comes from the zero eigenvalues. The second term comes from the negative 
eigenvalue. The third term is the determinant of the continuous positive eigenvalue spec- 
trum where det' means that the zero and negative eigenvalues are to be omitted. With 
these eigenvalues omitted the perturbed operator in the numerator has five less eigenval- 
ues in the spectrum because five of the eigenvalues of the unperturbed operator in the 
denominator have become part of the discrete spectrum of the perturbed operator. As- 
suming these eigenvalues have originated from the bottom of the unperturbed continuous 
spectrum we divide through the third term by a factor of uj"^ for each omitted eigenvalue. 
The contributions from the zero and negative eigenvalues are normalized with a factor of 
uj keeping each of the three terms dimensionless. 



4.1 Positive Eigenvalues 

The contribution of the positive eigenvalues requires the evaluation of the determinant 

dei'[-d^d^^ + f/"(06)] 
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However, since the {lo'"^)^ i^ ^^^ denominator corresponds to a set of measure zero in 
the continuous eigenvalue spectrum we can omit it and the notation det' which indicates 
omission of a discrete set of eigenvalues. 

This determinant ratio will turn out to be infinite and to obtain a finite answer we 
must divide by an infinite factor. The determinant in the numerator can be expanded in 
the following way: 
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It should be noted that the determinant is equal to the partition function of the self 
interacting massive scalar particle and that the second factor in the last line is expanded up 
to one loop contributions with three interactions with the effective external potential(see 
Fig. 1^). Tracing over a Cartesian basis we can see that the one and two interaction 
contributions are divergent but the three interaction contribution is finite: 






Figure 3: Expansion of the partition function. 
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where Vpert{k) is the Fourier transform of the Vperti'f')- We will only get a finite answer 
if we divide through by the infinite factors. 

The actual calculation is most easily done using hyperspherical coordinates in four 
dimensions, (r, 6, (p, ip), and expanding the eigenfunctions xi^^ ^) 0) V") i^^ terms of the 4D 
hyperspherical harmonics: 
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which are discussed in the appendix. 

In this situation the Laplacian, when operating on each term of 
number 'n', becomes: 
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There are 2n(n+2) + l such terms corresponding to different values of '1' and 'm' but with 
the same eigenvalue. Therefore: 
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First consider the denominator of (pSj), det[— P„ + u;^]. In order to calculate this 
determinant we need to solve the eigenvalue equation: 

[P„ - cj2 + \]u{r) = 0. (33) 

The solutions to this differential equation that are well behaved at r = and r = oo are: 
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Figure 4: The "potential" U"Mr)). 



where ux = \/\ — u"^. jn are the 4D analogs of the spherical Bessel functions and J„+i(x) 
are Bessel functions of the first kind. These solutions become purely oscillatory as r ^ oo. 
Notice that this solution is only well defined for A > u;^ and indeed there are no solutions 
for values of A < a;^ which are well behaved at r = and r = oo. Therefore the continuous 
spectrum of eigenvalues can be written as A 



u"^ + ul and: 
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The numerator of (pSj) involves the "potential" U"{(j)b{r)). For the symmetrized ef- 
fective potential (|T^ the "potential" for this problem is approximately constant except 
for a small perturbation in a small region near the bubble wall (see Fig.|^). The constant 
potential has been analyzed above. The solutions with or without the perturbations are 
identical for r -C i? and almost identical for r 3> -R. The perturbed solution differs from 
the unperturbed solution in this latter region at most by normalization and a phase shift, 
CbxT = uoxT + 5{ujx)- In this situation, for each value of n, we can obtain the ratio of 
determinants for a discrete spectrum by p6| : 
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which becomes in the continuum: 
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This formula gives the determinant ratio for a particular value of 'n' if we know the 
phase shifts, S{ujx)- In order to calculate the phase shifts, however, we must make further 
approximat ions . 

Consider both the perturbed and the unperturbed equations: 
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Multiply (^) by v{r) and (|39D by u{r), subtract, and integrate from to oo to obtain: 
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Using the fact that both solutions vanish at r = and using the asymptotic form of the 
Bessel functions we obtain an exact formula for the phase shift, S{ujx)'- 
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Of course, since the perturbed differential equation is extremely difficult to solve, we use 
perturbation theory to obtain: 



S{uJx) ~ 5{uJx)o = — 



OO yr rOO 

, yperti'^M^)^dr = - Vpert W^'^n+i(wAr)^dr, 

LOX JO Z Jo ^ 



(42) 



for small phase shifts. Using this result for the phase shift in Equation (|37D, we obtain 
the ratio of determinants: 



•det[I?„ + t/"(0,)] 



det[-Vn + uj^] 



exp 



oo /-oo UJ\T 

duj\dr^^ ^^pert('^)^n+i(^Ar)^ 



JO 



u^ + uji 



(43) 



for each value of n. The complete determinant ratio then becomes 



det[-9^9^ + cj2] 



exp^(2n(n+2)+l) 



n=0 



oo /-oo ^^^ 

io ^"^^^^:;^T^^p-t' 



r)Jn+i{uJxrf 



(44)- 

for an arbitrary perturbation potential. The next step in our analysis is to take into 
account the specific potential in our problem. 

It should be noted that this approximation gives the same result as the exact answer 



expanded to one loop (see Eq.^ and ^^ where the trace is taken instead over |rn/m), 
\knlm) bases. 

This formula is suitable for numerical evaluation, but in order to obtain an analytical 
answer we must approximate the perturbation in the potential U" by: 



K 



pert' 



10(r-M)2 g^ 







10 



^ > r > R 
10 ^ ^ '' 

otherwise 



(45) 
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Figure 5: The perturbation Vpert{r) is the dotted hne and the approximate perturbation 
Vperti''^) is the sohd hne. 



where P = (t/"(0(O)) — t/"(0(i?)))(see Fig. |^). With this approximate potential we find 
the phase shift via (^21): 



Soiujx) = -200-^ 



rR {Jn+i{ruJx) - Jn{ruJx)Jn+2{ruJx) 



121 
200^ 

^ (n + 4)!(n+l)! 



2, ,2 



2{n + 2)(r2 + 4) 1F2 n + -, ra + 4, 2n + 3; -r^w 



-r^tu^ 1F2 ( n + -, n + 5, 2ra + 4; -r^w^ 



112-2-2nr3i?(rcuA)2+2" .35 7 

; ^T77 ,.T„ 2-^3 ri H — ,n -\ — , n + 2, n H — , 2?t, + 3; — r 0;,, 

10 (n + |)[(n + l)!]2 ^ n 2' 2' ' 2' ' ^ 



11/ioK 



.(46) 



i? 



Using this result we can verify that S{uJx) -C 1 which means that the approximations 
of ( ^2] ) are justified. The phase shift as a function of ux is shown for a few values of n in 
Fig.(^). However, (^61) is very difficult to work with so we make another approximation. 
For ux large enough we can use the asymptotic approximation for the Bessel function in 



Jn+lix) = \ <^ Pn+lix) COS 

7TX 



3\ TT 

X — [ n -\ — — 
2/2 



3\ vr 

X — [ n -\ — — 

2/2 



where: 



P„+i(a;) = 1- 



-Q„+i(x)sin 
(/i-l)(/x-9) , (^-l)(^-9)(/x-25)(/i-49) 



2! r8x)2 



4! fSxl^ 



, (47) 
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Figure 6: The phase shift 6q{ujx) shown as a function of ux for various values of n. Notice 
the changing vertical scale. 



Qn+i{x) 



(/i-1) (/,-l)(;,-9)(/i-25) 



/i 



In this way we find: 



1! {8x) 
4(n + l)l 



3! (8x) 



+ ■■ 



5o(^a) 



2RP 



+ 



-l)"/?COs(2i?CJA) 



TT 



UJ\ 



+ Oil/iv, 



(49) 
(50) 

(51) 



Of course this estimate for 6o{ujx) has an infrared cutoff in 0;^ below which it is not 
a good estimate and this will be true no matter how many terms in the approximation 
of the Bessel function we keep. Therefore as a first approximation we keep only the first 
term and cut off the ux integral at the position of the first peak in 6q{ujx). Using this 
approximation for 6q{ujx) we obtain our first estimate of the ratio of determinants for each 
value of the quantum number 'n' given by: 



det[-Vn + U"{<f).)] 



exp 



2R(3 f u 

arctan 



SOttcj 



uj{n) ^ 



(52) 



where uj{n) ^ 33.16ra + 72 = an + 6 is the infrared cutoff given by the location of the 
first peak in 5o{ujx)- Using the result (B2|) we can calculate our first approximation for the 



complete determinant ratio. For large values of n the terms of the sum in (^4]) approach: 



uo 



4RP 2 

n 

30ttuj \an 



(53) 
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and therefore the sum is infinite. Our first estimate for the complete determinant ratio 

based on (|^) is infinite. This was to be expected based on comments at the beginning of 

this section and, since this approximation exactly coincides with the tadpole diagram, we 

subtract this term from the exponent as our normalization prescription for the complete 

determinant ratio. Note that we have used the approximation of ([51]) to calculate the 

infinite factor but the equality with the tadpole diagram holds before this approximation 

was made. 

In order to find a finite result we must adjust our earlier approximation for the phase 

shift: 

1 3 

6{uJx) = arcsin((5o(^A)) = ^oK) + tt-^oK)^ + ttt^oI^a)^ + ■ ■ ■ (54) 



which would add correction terms to the exponent in (|5 
to a factor of the form: 

LUxSoiuJxY 



exp 



IT 



YaL "^"^''co^ + col 



The second term in (|5^ ) leads 

(55) 



for each value of n. This term is also divergent but does not exactly coincide with the 
next term in the expansion (|26|) . However, the divergent contribution in each must be 
the same. 

Therefore dividing through by the previously obtained infinite factor [ie. the right 
hand side of (0)) leads to the renormalized value: 



det[-Vn + U"i(l)i>)] 
det[-Vn + U"{<j)^)] 



exp 



TT 



(2R(3\ 



ijj 



24 \307rcj/ \^{n) 



arctan 




(56) 



For large values of n the terms of the sum in ( ^41) approach: 



7r2 / 2i?/3 \ 



24 



Uottcj J 



(2n^ + 4n + 1) 



LO 



Za^rfi 



(57) 



which also leads to a divergent sum and therefore an infinite value for the determinant 
ratio. The divergence here must be contained in the divergence in the two interaction 
term of the expansion (p6|) , as are other logarithmic divergences obtained from keeping 
more terms in the approximation of the Bessel function (P7|). However, subtracting the 
two interaction term will almost certainly leave a finite contribution at the next order in 
n. We will assume this is the case but since we are doing an approximate calculation we 
will not calculate the finite contribution from the two interaction term or any correction 
terms to our approximation since they will have the same physical dependence as the 
approximation we will give. We should further note that correction terms are probably 
not calculable analytically and are not likely to be a problem for our results. They are 
obtained from integrating oscillatory functions (see second term in (0)) over io\ in (|55|) 
and should not contribute very much. The renormalized approximation to the complete 
determinant ratio (|^) that we obtain from (|56D is given by: 



exp 



24 VaOvrcuj 



^(2n(n + 2) + l) 



ix) 



.n=0 



uj\n] 



arctan 




E2 

n=l 



UJ 



n 



?)d^n^ 



(5^ 



15 



-4(Sli^«("+^'+^'{ 



100 



uo 




E { (2n(n + 2) + l)£^^\ + ^C(2) + £,C(3) 



n=l 



Sa'^n^ 






Evaluating gives: 



exp 



7r2 2R13\ 

-5300— 

24 VSOttcj/ 



exp [-2 X 10 



-5l 



(59) 



(60) 



Because of the small value of the exponent we can see that the exponential is extremely 
well approximated by: 



exp 



-53007r2 / 2R(3 



24 VSOttcu 



1- 



53007r2 / 2R[3 " 
24 VSChRu^ 



(61) 



which is our result for the determinant ratio (p5|). 

This is a nonperturbative calculation because Eq.(|37D is a nonperturbative resum- 
mation of the perturbation expansion ([26| ) of the determinant. While we make an ap- 
proximation through Eq.(^) this expansion is nonperturbative since the terms of the 
expansions do not coincide. Our approximate calculation should therefore have contribu- 
tions from all perturbative diagrams. This nonperturbative resummation of diagrams in 
the determinantal prefactor is very similar in spirit to that of |2^ . 

This result constitutes the contribution of the positive eigenvalues to the quantum 
corrections to the decay rate in the zero temperature theory. In the next section we 
consider the zero and negative eigenvalues. 



4.2 Zero and Negative Eigenvalues 



The zero eigenvalues contribute J B/2ti per collective coordinate P]. The action of the 

4D bubble is independent of the center of the bubble which means there are 4 collective 

coordinates leading to the first factor in Eq. 

The eigenfunctions of zero eigenvalue are: 



Xopi\x., y, z, t) 



d 



<Pb{r) 



dr d 



dx'^ dr 



Mr) 



xf" d 
r dr 



<Pb{r) 



where r = a/x^ + y'^ + z'^ + t''^ and: 



t = 


r Cosiij) 


~ '^100, 


z = 


r Cos{e)Sin{ilj) 


~ Yno, 


y = 


r Cos{(j))Sin{6)Sin{il)) 


~ Yin + Yii-i, 


y = 


r Sin{(j)) Sin{6) Sin{ilj) 


~ Yiii — Yii-i. 



(62) 



(63) 
(64) 
(65) 
(66) 

Therefore these eigenfunctions all correspond to ra = 1 and since there are no radial nodes 
we can be sure that there are no negative eigenvalues with n 7^ 0. 
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In |2B|, Coleman argued that there is only a single negative eigenvalue for an 0(4) 
invariant bounce. We assume this^ and determine its value. We use the method of 0] 
with a slight modification. As was argued by Coleman, the only possible eigenfunctions of 
negative eigenvalue are those that are bound to the bubble wall. For such eigenfunctions 
we can approximate the centrifugal potential in ( |3TD by a constant determined by its value 
at the bubble wall {r = R): 



X 



pn 



4n2 + 8n + 3 



(67) 



where Xp is a number independent of n. We know that for n = 1 the lowest eigenvalue is 
zero: 



A, 



01 



Ao + 



15 

IR2 



0. 



Therefore we can obtain the lowest eigenvalue for n = 0: 



A_ = A 



00 



Ao + 



15 



4i?2 



4i?2 4i?2 



3 



{61 



(69) 



This value is different from Coleman's but only by a factor of 2. We cannot explain this 
discrepancy but can only stand by our calculation. 

4.3 Decay Rate for Zero Temperature 



Therefore combining the zero temperature semiclassical result (jT^) with the quantum 
corrections obtained in the last section we estimate the decay rate per unit volume of the 
false vacuum in the zero temperature theory to be: 



r/v 



exp 



■^'' K (i; 



UJ 



v^ 



det'[-d^d'^ + U"{<j)b)] 



(io-^)^ det[-^f,^^' + 00^] 



exp 



277rVn /277rVcu2^ 



Rio 

'2{AEy) [47i{AEf) 71 

1.55 X 10~^^MeV^ ^ 3 x lO^^fm^^s 



53007r2 (2RP\- 
24 \307vuj) 



-1/2 



(70) 

(71) 
(72) 



It should be noted that the quantum corrections are negligible so the decay rate is basi- 
cally determined by the semiclassical result. We believe that this observation would be 
unchanged by the inclusion of corrections that we have ignored in this calculation. While 
the quantum corrections did not turn out be significant in this case, the fact that they are 
not is relevant to the baryogenesis mechanism described in 0. As well, the techniques 
applied to the problem may prove useful in other calculations of this type. 

In the next section we perform the same calculation in the high temperature limit. 



^This is most likely a good assumption but not proven due to the cusp in our potential. For more 
details see §^,§^ and IM 
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5 Quantum Corrections for High Temperature 



The quantum corrections to the decay rate at high temperature correspond top3| 



A = T\u^ d^x 



2nT J 



det'[-did^ + U"{<P,)] 



{uj-^)^ det[-did^ + U"{(j)-)] 



-1/2 



(73) 



which again factors into three parts corresponding to zero, negative and positive eigen- 
values respectively. 

5.1 Positive Eigenvalues 

The calculation in three dimensions is extremely similar to the four dimensional calcula- 
tion presented in the Sect. 4. The expansion of the determinant ratio is exactly the same 
as in 4D (|26|) . Tracing over Cartesian bases we can see that the one loop term is divergent 
while the two loop term is finite: 




d'x d^'k ^^^^'^""^ 



fc2+a;2' 



j3;, j3„ ypertyk) Vpert{ k) 



d k dp 



(A;2 + a;2)((fc + p)2 + o;2)' 



(74) 
(75) 



The calculation is done using spherical coordinates and the eigenfunctions are ex- 
panded in terms of the spherical harmonics: 



oo n \l\ / \ 

n=Ol=-nm=-\l\ ^ 



(76) 



For each value of / the Laplacian becomes: 
Therefore: 



d'^ l{l + 1) 
(ir2 r2 



r = —Vir. 
r 



det[-did' + U"] = n(det[-P, + U"]) 

1=0 



f/l\2l+l 



(77) 



(7J 



The solutions to this differential equation that are well behaved at r = and r = oo are: 

u{r) = ujxrji{ujxr) = ^J^-Ji+i/2iuJxr), (79) 

where ji are the usual spherical Bessel functions and tux = VA — cu"^. We obtain: 



1 f°° 
6{ux) ^ — / K)ert' 
ujx Jo '■ 



TC 



r)u(r) dr = — 



2 Jo 



ypeitir)rJi+i/2{uJxrydr, 



^0) 
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thus givingf^l^: 

det[-did'^ + U"i(l)b)] 



det[-did' + 00^] 



exp^(2/ + l) 






JO 



UJ^+LOI 



1+1/2 



(uJxry 



(81) 

As in the zero temperature case we can find the exact expression for the phase shift with 
the approximate potential but it is too difficult to work with. Instead we approximate 
the Bessel functions as in {^7\ ) where {n + 1) —>■ {I + 1/2). In this way we find: 



6{ujx) ~ So{uJx) 



'"'' +0(1/-;) 



SOnux 



(82) 



Our first estimate of the ratio of determinants for each value of the quantum number '1' 
is given by: 



det[-A + f/"(0b)] 



exp 



SOvru; 



arctan 



U) 



u!{n) 



(83) 



det[-P, + f/"(0„)] 

where ui{l) ~ 49.5/ + 82 = d + d is the infrared cutoff given by the location of the first 
peak in 6o{ujx) as a function of ujx- For large values of / the terms of the sum in 

approach: 

4R/9 /,,,\ 



ARI3 fuo 



and the sum is infinite and therefore the complete determinant ratio is also infinite. Again 
this was to be expected and, since this approximation exactly coincides with the one loop 
term, our renormalization prescription is to remove the one loop term from the exponent 

ofdiD. 

Again we must adjust our earlier approximation for the phase shift: 



1 3 

5{ujx) = aTcsm{5o{uJx)) = 5o{ux) + -5o{uJxY + -jT^^oi^^xf + 

o 4U 



^5) 



and the correction is implemented as in (55). Evaluating the contribution to the deter- 
minant ratio leads to: 



det[~Vi + U"{<P,)] 
det[-Vi + U"{<f)^)] 



exp 



n^ / 2Rp y / u 
24 VaOvrcj/ [Oil) 



— arctan 



u 



Mi)^ 



(86) 



For large values of / the terms of the sum in (^) approach: 



IT 



( 2RI3 



24 VSOtto;; ^^^^^^(sc^pj- 



^7) 



^°Note that we have again dropped the factor of {oj^'^Y ^-^d the "det"' notation as the omitted eigen- 
values correspond to a set of measure zero. 

^^The summand in this expression appears in [|l6[ but is only used for the high energy modes. The 
non-divergent term in what follows does not appear in |l6[ . 
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The sum is finite so we do not need to further normahze and our result for the complete 
determinant ratio (0) is: 



exp 



TT^ / 2Rf3 ' 

24 \30nuj^ 

100 c 

-E (2^ + 1) 
1=1 I 



100 ( 

E(2^ + i) 

1=0 I 



— arctan 



3 ^ 2uj^ 



m. 



3^hi^^(^) + l^^(^) 



Evaluating gives: 



exp 



24 \307rujJ 



exp[3 X 10 



-7i 



^9) 



Because of the small value of the exponent we can see that the exponential is extremely 
well approximated by: 



exp 



827r^ ( 2R(3 \ 

24 \m^) 



1 + 



827r2 



( 2Rp \^ 



24 \ SOvrw / 



(90) 



This result constitutes the contribution of the positive eigenvalues to the quantum cor- 
rections to the decay rate in the high temperature theory. In the next section we consider 
the zero and negative eigenvalues. 



5.2 Zero and Negative Eigenvalues 

In the 0(3) invariant bubble there are 3 collective coordinates leading to the first term in 
([73|). The eigenfunctions of zero eigenvalue are: 



where r = \/x^ + y'^ + z'^ and: 



d 



Mr) 



dr d 



x^ d 



dx^^ dr r dr 



(91) 



z= r Cos{9) 


~ ^10, 


y= r Cos{(p)Sin{9) 


~ Yn + Yi_i 


y= r Sin{(p)Sin{e) 


~ Fii - Fi_i 



(92) 
(93) 
(94) 

Therefore these eigenfunctions all correspond to Z = 1 and since there are no radial nodes 
we can be sure that there are no negative eigenvalues with / 7^ 0. 

We will assume, as we did in the previous section, that there is only a single negative 
eigenvalue and concentrate on obtaining an approximation. The only possible eigen- 
functions of negative eigenvalue are those that are bound to the bubble wall. For such 
eigenfunctions we can approximate the centrifugal potential in ([77|) by a constant: 



A 



pn 



Ap + 



/(/ + 1) 

i?2 ■ 



(95) 
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where Xp is a number independent of /. We know that for / = 1 the lowest eigenvalue is 
zero: 

Aoi = Ao + |j = 0. (96) 

Therefore we can obtain the lowest eigenvalue for / = 0: 

A_ = Aoo = Ao = -^. (97) 

5.3 Decay Rate for High Temperature 

Therefore combining the high temperature semiclassical result (^2D with the quantum 
corrections obtained in the last section we estimate the decay rate per unit volume of the 
false vacuum in the high temperature theory to be: 






827r2 f2Rl3\ 



-1/2 

• (98) 



We cannot obtain a numerical estimate for the decay rate at this time since the temper- 
ature dependence of the semiclassical result is not currently known. This result could be 
important to the study of false vacuum states of this type in heavy ion collisions. 

6 Conclusion 

We have obtained an estimate for the decay rate per unit volume of T /V ^ 3 x lO^^fm^ s^^ 
for a false vacuum in an effective Lagrangian approach to QCD for zero temperature 
theory. We have obtained an expression (p8|) for the decay rate per unit volume in the high 
temperature theory, which is the best we can do without knowledge of the temperature 
dependence of the effective potential. These are nonperturbative calculations in that they 
contain contributions from all orders of perturbation theory. 

The value of the decay rate for zero temperature shows that if the universe had started 
out in a false vacuum state then it would have decayed long ago into the true vacuum 
state. This is not new or interesting. The interesting thing concerns the possibility that 
these vacuum bubbles would have significant effects on the evolution of the early universe 
at around the time of the QCD phase transition. 

One such effect related to baryogenesis was described in 0. The first approximation 
to the decay rate we have calculated shows that bubbles of false vacuum are far too 
unstable for this simplified baryogenesis mechanism. The effect of fermions or intrinsic 
heavy degrees of freedom could go a long way to stabilizing the bubbles because they 
would make the barrier separating the false vacuum from the true vacuum much higher, 
thus increasing the stability of the false vacuum. The precise calculation is outside the 
scope of the present work and we can only say that, so far, this mechanism for baryogenesis 
at the QCD scale has not been proved viable. 

The result does not, however, rule out the possibility of nontrivial effects on the 
evolution of the early universe. As well, the techniques developed may be useful in the 
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study of arbitrary metastable vacua which seem to be a general feature of gauge theories 
in the strong couphng regime. 

The high temperature expression for the decay rate will be relevant to the possibility 
of observing CP odd bubbles at RHIC through signatures as described in |T^. As was 



mentioned previously the false vacuum described in [0] assumed the large Nc limit. 
However, the same effects would manifest themselves in our case which arises from an 
effective potential valid for arbitrary Nc 

In the future it might be useful to verify numerically our contention that the correc- 
tions neglected in the approximation and finite contributions from the renormalization 
prescription do not significantly alter the results. As well, one could do numerical cal- 
culations using the exact numerical domain wall solution in the case of non-degenerate 
vacuum states. It is possible that the determinant could also be estimated numerically in 
the non-degenerate case. 



It would also be interesting to repeat this calculation using the methods of []TB| [|T^ 



Tql |T9|| and compare with the present results. 



A Hyperspherical Harmonics in Four Dimensions 

Hyperspherical coordinates in 4D dimensions are related to the Cartesian coordinates by: 

X = r Sin[ilj]Sin[9]Sin[(f)], (99) 

y = r Sin[ilj]Sin[e]Cos[(j)], (100) 

z = r Sin['^]Cos[e], (101) 

w = rCos[ij]. (102) 

The Laplacian in 4D in hyperspherical coordinates is: 
^dr (r'dr) + ^^d^ (sin^ ^d^) + ^ .1 de (sin Ode) + , . ^V ^ A ' ^^0^) 

Assuming separable solutions and treating 9 and coordinates in exactly the same way 
as in three dimensions we obtain the differential equation: 

^"(V^) + 2 cot(^)^'(V^) - /(/ + 1) csc2(^)^(V^) = K^i^l)). (104) 

With the substitution u = Cos[ip] this becomes: 

(1 - u^)U"{u) - 3uU\u) - ^-^^^U{u) = BU{u). (105) 

\i B = n{n + 2) and /(/ -|- 1) = /'(/' + 1) this can be identified as the differential equation 
satisfied by the associated type II Chebyshev functions. These can be obtained from the 
type II Chebyshev differential equation in exactly the same way as associated Legendre 
functions are obtained from the Legendre differential equation. 

As an aside notice that the type II Chebyshev equation is a special case of the Geigen- 
bauer (Ultraspherical) equation: 

(1 - x')^C\^\x) - (2a + l)x^C\^\x) - n{n + 2a)C^^\x) = 0, (106) 
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with a = 1. The Legendre polynomial equation corresponds to a = 1/2. All other 
hyperspherical coordinates will lead to associated Geigenbauer equations with integer or 
half integer a. 

The hyperspherical harmonics in four dimensions are given by: 

V (ft rh ,IA - A(n I m\S '^lrn{0,(p) [/n^(cOs(^)) 0<l<n 

where Yim{0, 0) are the usual 3D spherical harmonics and Uj are associated Chebyshev 
type II functions defined by: 

UJ{X) = (l-x^y/'-^Unix) 



(l_^2y/2.(-ir(« + l)V^c?' 



2"+i(n + l/2)! dx^ 
for / > and by: 






,(108) 



I ( l)"(n+l)v/^ _ 2y/2 ^"^'^^ Tq _ 3,2Nn+l/2l QggN 

for / < —2. These hyperspherical harmonics form a complete orthogonal basis for the 
functions of the angular variables in four dimensions. 

E E E Y:,„Xr,eA,mnUr,e',<j>',^') = ^-^^i^^^-^^^6{<j>-<j>'). (no) 

~Ui , sm -w smb* 

n=Ul=—nm=—l t 

oo riv ("K r2n 

/ / / Y:i^{r,e,(j),^)Yr,n'm'{r,e,(j),^) = Snn'Si,6u'. (Ill) 

Jo JO JO 
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